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Summary: 


^The  general  equations  of  motion  for  a  multicomponent  chemically 
reacting  flow  system  are  introduced.  The  concept  of  influscid  flow 
is  presented  and  the  fundamental  boundary  layer  equations  appropriate 
for  the  analysis  of  the  interaction  between  a  flowing  multicomponent 
chemically  reacting  gas  and  a  surface  are  formulated  for  a  typical 
coordinate  system.  (  )  ^ 
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1.  INTRODUCTION 


In  investigating  the  interaction  between  a  surface  and  its  environment, 
one  is  concerned  with  an  auialysis  ol  fluid  dynamic  fields  involving  the 
simultaneous  transport  of  mass,  momentum  and  energy.  Such  problems 
are  described  by  the  equations  of  change  which  are  a  coupled  system  of 
nonlinear  partial  differential  equations.  The  equations  of  change  comprise 
the  conservation  laws  for  the  fluid  system  and  consist  of  the  conservation 
of  chemical  species^  the  momentum  equation  and  energy  equation.  The 
global  form  of  the  conservation  of  mass,  or  continuity  equation  can  be 
obtained  directly  by  summing  over  all  of  the  equations  for  the  conservation 
of  chemical  species. 

The  presence  of  chemical  reactions  in  multicomponent  gaseous  flow 
systems  greatly  complicates  the  analytical  formulation  and  solution  of  the 
flow  field  problem.  The  complexity  of  the  problem  and  the  relatively 
small  effort  expended  in  this  area,  as  opposed  to  the  voluminous  literature 
and  effort  dealing  with  classical  nonreacting  boundary  layers  and  inviscid 
flows,  have  conspired  to  limit  the  current  state  cf  our  knowledge  concerning 
this  technically  important  class  of  problems.  It  is  therefore  of  considerable 
interest  to  consider  the  general  equations  of  change  appropriate  to  a 
multicomponent  chemically  reacting  flow  system,  and  from  these  to  derive 
the  boundary  layer  equations  which  describe  the  phenomena  in  the 
immediate  vicinity  of  the  surface  of  a  body  immersed  in  such  a  fluid. 
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An  interesting  point  is  that  in  the  classical  treatment  of  flow  fields, 
the  problem  is  often  split  into  the  solution  of  two  problems,  inviscid  and 
viscous.  Ail  analogous  treatment  for  reacting  fields  would  be  to  define  an 
”influ8cid**  flow  field  as  that  region  where  gradients  of  velocity,  tempera¬ 
ture,  pressure  and  species  concentration  are  so  small  that  transport  pro¬ 
cesses  related  to  viscosity,  thermal  conductivity,  thermal  diffusion,  pressure 
and  concentration  diffusion  may  be  conveniently  neglected.  The  **  fluscid" 
region  is  then  the  boundary  layer  adjacent  to  the  body,  where  the  afore¬ 
mentioned  gradients  are  so  large  that  coupled  transport  processes  constitute 
the  essential  physical  phenomena.  It  is  noted  that  in  the  case  of  a 
reacting  flow,  in  addition  to  gradients  of  velocity  and  temperature,  species 
concentration  gradients  also  appear  near  the  body  surface,  due  to  the  presence 
of  dissociation,  combustion  and  other  forms  of  chemical  reaction.  Conse¬ 
quently,  one  must  regard  the  boundary  layer  as  a  multidimensional  region. 

The  influscid  and  fluscid  fields  can  then  be  solved  separately,  and 
matching  is  satisfied  by  taking  the  ** inner"  boundary  conditions  for  the 
influscid  field  as  the  conditions  as  "infinity"  for  the  boundary  layer.  Thus, 
in  simplifying  the  analysis,  by  reducing  it  to  the  solution  of  two  problems,  a 
complication  is  introduced  in  the  form  of  interaction  between  the  influscid 
and  fluscid  fields.  In  the  case  of  a  classical  non-reacting  gas  flow,  when 
the  inviscid  flow  field  is  irrotational,  the  interaction  appears  primarily  as  a 
displacement  effect  at  the  wall,  in  the  form  of  a  thicker  body  or  equivalent 


boundary  layer  displacement  thickness.  However,  when  the  inviscid  flow  is 
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rotational,  as  for  example,  when  there  are  entropy  gradients  in  the  field  due 
to  curved  shocks,  then  in  addition  to  the  displacement  effect,  there  is  an 
interaction  between  the  vorticity  of  the  external  field  and  the  vortical 
boundary  layer.  For  the  case  of  weak  interaction,  it  appears  that  this 
effect  may  be  included  in  the  boundary  layer  analysis  by  modifying  the 
boundary  conditions  for  the  velocity  gradient  at  infinity^. 

In  a  reacting  floW|  another  type  of  interaction  comes  into  play  in  which 
the  concentration  gradients  must  be  made  to  match  at  the  outer  junction  between 
the  fluscid  and  influscid  solutions,.  And  further,  one  cannot  now  take  the 
surface  conditions  independently  of  the  surface  mass  and  energy  balance, 

Int'  itively,  one  would  expect  that  all  of  these  co  .pling  effects  could  be 
treated  by  means  of  an  iterative  approach. 

In  considering  the  flow  of  a  reacting  gas  over  a  surface,  it  is  not 
correct  in  general,  to  treat  the  environment  as  a  single  gas  unless  at  each 
point  the  diffusion  flux  set  up  by  concentration  and  thermsJL  gradients,  and 
the  convective  flux  set  up  by  pressure  gradients ^  are  exactly  counterbalanced 
by  the  influx  of  species  by  coxivectioxi  and  the  production  of  species  by 
chemical  reactions,  so  that  a  stationary  equilibrium  composition  can  be 
maintained  in  the  flow  field.  Hence,  if  exact  solutions  are  desired,  when 
chemical  components  are  free  to  react,  expecially  when  a  new  species  is 
introduced  into  the  boundary  layer  by  means  of  mass  transfer  from  the 
surface,  due  to  vaporisation,  sublimation,  heterogeneous  reaction,  or 
direct  fluid  injection,  it  appears  necessary  to  evaluate  the  transport 
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coefficients  and  thermodynamic  properties  including  viscosity,  thermal 
conductivity,  concentration  diffusion,  thermal  diffusion,  specific  heat 
and  enthalpy  at  each  point  in  the  boundary  layer  as  a  function  of  the 
local  composition,  pressure  and  temperature.  Moreover,  it  is  essential 
that  the  driving  forces  for  each  species  be  determined  by  means  of  the 
individual  equations  for  the  conservation  of  species,  so  that  the  extent  of 
surface  reactions  can  be  assessed. 
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2.  SYMBOLS; 


CL  acoustic  velocity 


Bi-  — 

^  n 


referred  mole  fraction  of  i  species 


Cl  mass  fraction  of  species  i 

c,  ..  Pi  .  = 

■  7/1 


(2  0,-0 


^  » 

Cp  specific  heat  at  constant  pressure^  Cv  specific  heat  at  constant 
volume 


driving  force 


diffusion  coefficient»  (multicomponent  system) 


diffusion  coefficient,  (binary  system) 


thermal  diffusion  coefficient 


thermodynamic  internal  energy  per  unit  mass,  including  chemical 


E  stagnation  internal  energy^  E  =  C  +  —  XT  ^ 

F:  external  force  acting  on  a  unit  mass  of  species  i 

h  static  enthalpy ,  h  *  S  +  _£ 

P 

static  enthalpy  of  the  species,  including  enthalpy  of  formation 

*■  Pi 
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A\n 


.th 


r  enthalpy  of  formation  of  i  species 


H  stagnation  enthalpy ,  H  -  x  V” 


1 


unit  tensor 


^  i  *  relative  mass  flux  of  species  i 

absolute  mass  flux  of  species  i 
relative  molar  flux  of  species  i 


A 


coefficient  of  thermal  conductivity 


BoltsmannU  constant 


/€  linear  dimensicni  cm* 

molecular  weight  of  species  i 

molecular  weight  of  the  (;as  mixture 

e 

Mi  net  molar  production  of  species  i  per  unit  volume  by  chemical 
reaction 


moles  of  i  per  unit  volume 

total  number  of  moles  per  unit  >  olume 

number  of  molecules  of  i  per  unit  volume 
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total  number  of  species 


static  pressure 


energy  flux  vector 


universal  gas  constant 


radius  of  curvature 


radius  I  measured  from  the  centerline  to  the  surface 


time 


temperature 


X  component  of  velocity 


y  component  of  velocity 


absolute  flow  velocity  of  species  i 

-weighted  average  velocity  of 
"Vj  •  Vj,  -  diffusion  velocity  of  species  i 


the  fluid  mixture 


net  mass  rate  of  production  of  species  i  per  unit  volume  by 
chemical  reaction  where  S  •  o 

I  ^ 


mole  fraction  of  species  i  , 


Xi «  !1l 


i  ZXi 

i 


-  8  - 


coordinate  system 
symmetric  rate  of  strain  tensor 
viscous  stress  tensor 
ordinary  viscosity  coefficient 

a  quantity  related  to  the  second  viscosity  coefficient 

partial  density  of  species  i, 

density  of  the  fluid  mixture,  ^  «  Y\^7lfl 

pressure  tensor 

dissipation  function 


Subscripts 

i  i^  chemical  species 

•b  total 


Dimensionless  Groups 


Prandtl  Number 


Lewis  Number 


Thermal  Lewis  Number 


3.  BASIC  RELATIONSHIPS 


Referring  to  the  list  of  symbols ,  we  observe  that  the  total  number 
of  moles  in  a  unit  volume  is  obtained  by  summing  over  all  the  contributions 
of  each  species  I  i.e. 

2  Y\i  (1) 

1 

The  density  of  the  fluid  is  given  by  summation  over  all  the  partial 
densities  of  all  the  species. 

•  S  Pi  ‘  “  2./3Ci  (2) 

The  absolute  molar  flux  of  species  i  with  respect  to  fixed  spatial 
coordinates  is  given  by: 

niVi 

hence  the  mass  flux  of  species  i  with  respect  to  fixed  spatial  coordinates  is 
given  by; 

3i  •  ■  fiXi  (3) 

We  may  thue  obtain  the  maee-weighted-average  velocity  of  the  fluid 
by  summing  over  Ji 
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Thi«  velocity  V  ie  the  mecroecopically  observed  etreem  velocity*  i.e. 


jr  ’  ^  2  li  •  ^  E  -  i  I •  2ciVi  (5) 

The  diHation  velocity  ^ i  oi  epeciee  i  ie  then  defined  ae  the 

difference  between  the  absolute  velocity  of  species  i  and  the  obsefved  stream 


velocityo 


Vi  *  Vi-v- 


It  is  then  possible  to  define  the  diffusion  flux  vector 

ji  •  /“lYi  p) 

Note  that  the  mass  weighted  average  of  the  diffusion  velocity  is  zero. 
This  may  be  shown  quite  easily,  by  summing  equation  (6)  over  all 
and  then  introducing  equations  (2)  and  (4). 

pS/’iYi  •  j;  IftYi- 


Consequently  we  may  write: 


Zji  -  ZpiY-0 


(8) 
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4.  CONSERVATION  OF  SPECIES 

The  conservation  of  species  i  may  be  written: 

^  +  v.(ni-iri’)  =  (9) 

or  equivalently 

^  +  V-Ii  »  Wi.  (10) 

where  the  chemical  source  function  represents  the  mass  rate  of 

production  of  species  i  by  chemical  reaction  and  may  be  determined 
only  with  great  difficulty  from  chemical  kinetics. 

Upon  eliminating  ^  t  we  obtain: 

+  V-  [nilv*Vi^]  .  Mi 

If  there  are  N  different  species  in  the  flow  fieldi  only  N*1  of 
these  equations  will  be  required  in  addition  to  the  global  continuity  equation, 
which  follows. 

Summing  equation  (10)  over  all  the  species,  we  find: 

=  zw, 

Since  mass  Is  neither  being  cresteci  nor  destroyed)  neglecting  relativistic 
effects, 

S  Wj  =o 
i 


(13) 
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and  the  former  becomes: 

^  v  -o  04) 

St 

the  global  oontlauity  equation.  Other  forms  may  also  be  derived. 

These  are  listed  below  for  convenience. 


±  t 

dt 


7-‘Vr  a  O 


P  7V 


dCi 


dt 


it  the  etreaming  derivative.  Undoubtedlyg  we  have  not  exhausted  all  the 
pottibilitiet  in  representing  the  conservation  of  mass. 

Returning  briefly  to  the  flux  with  respect  to  the  mass  average  velocity, 
we  note  that  the  flux  may  be  broken  down  into  several  parts,  i.e. 

ii'PiYi'/’i  4  y/"''  4  yi<f  >  4y,,  )  a^) 

Here  we  observe  that  the  mass  flux  ^  ^  contains  components  due  to  gradients 
in  concentration,  temperature,  pressure,  and  external  forces.  From  the 
thermodynamics  of  irreversible  processes,  if  the  situation  encountered  is 
not  too  far  removed  from  equilibrium,  the  flux  must  be  a  linear  function 
of  the  driving  forces,  and  the  net  flux  is  then  the  sum  of  the  ind'vidurl 
contributions.  Hirschfelder  etal^  gives: 
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\^ere 


Z  WiWjliijdj  -  vlnT 


(18) 


il  ■  7Xj  +(Xj  -Cj^vl^p 


-  £i  t  ^  F, 

p  l-  Wj 


2  *^K  F  K 


] 


For  the  cate  of  a  binary  mbcture  in  which  only  concentration  diffusion  and 
thermal  diffusion  are  aignificantt  equation  (IS)  reduces  to: 


ji  »  vXj  -  Di  vlnT 

and  it  it  then  not  difficult  to  show  that: 


(20) 


Yi  ■  -  tCf^  ^  vlnT  (^i) 

We  observe  that  no  simplifications  residt  upon  introduction  of  mass 
fraction  rather  than  mole  fraction*  when  there  are  more  than  two  chemical 


components. 


5.  THE  MOMENTUM  EQUATION 


The  momentum  equiUion  in  aerothermochemistry  is  identical  with 
that  used  in  aerodynamics  iof  gases  without  chemical  reactions. 

<“> 

«h.r.  is  Um  •jrt.riuU  fore*  .etlog  oa  «  unit  ru**  of  species  i  , 

MMi  *tT  is  the  preseure  tensor.  The  pressure  tensor  may  be  written; 

H  ■  -pi  +  III  (23) 

where  p  is  the  static  pressure  and  ^ is  the  viscous  stress 
tensor.  According  to  equation  (22)p  the  macroscopic  velocity  of  the  fluid 
undergoes  a  change  because  of  the  gradient  in  static  pressure,  internal 
shear  and  the  external  forces  which  act  on  the  various  chemical  species 
present.  We  may  write  for  ; 

lij  =  +  I>'^T  («) 


udisrs 


SiJ 


is  the  conventional  viscosity 

is  related  to  the  second  viscosity  coefficient 

is  the  unit  tensor 

is  the  symmetric  rate  of  strain  tensor  where  a  typical  term  is 
given  by: 


€ii  -  i  (1^  . 

'  SXj  ) 

•uch  that  ,  in  cartesian  coordinates. 
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If  it  is  assumed  that  Stokes*  postulate  holds  which  is  equivalent  to  assuming 
that  we  are  dealing  with  a  process  whose  characteristic  time  is  much 
omaller  than  the  viscous  characteristic  time*  then 

(26) 

and  there  follows: 

Tii'  2/*  (Si)  - 


(22) 


6.  THE  ENERGY  EQUATION 


-  16  - 


A  quite  general  form  of  the  energy  equation,  is  given  by: 

=  -v-g  < 

This  aquation  is  formally  analogous  to  the  first  law  of  thermodynamics, 


where  the  internal  energy 


includes  chemical  energy.  In  words, 


this  equation  etates  that  the  internal  energy  of  the  mixture  6  ,  following 

the  RMan  motion,  changes  due  to  the  combined  effects  of  the  energy  flux 
vector  Q  ,  the  work  of  the  pressure  tensor  ^  and  the  work  done 
by  the  external  forces  . 

The  energy  equation  may  be  rewritten  in  many  different  equivalent 
forms,  but  before  we  develop  some  of  these  other  forms,  it  will  be 
instructive  .to  examine  the  energy  flux  and  pressure  tensor  work  terms  a 
little  more  closely.  The  energy  flux  with  respect  to  the  mass  average 
velooity,  denoted  by  0  may  be  broken  down  into  its  component  parts,  i.c. , 


Cci^ 


which  repreient  contributions  due  to  temperature  gradients »  diffusion^ 
concentration  gradients  and  radiation*  If  the  effects  of  radiation  can  be 
separated  outg  the  remaining  termn  are^: 


Q  =  -KVT  +  XpiYjht  - 


llni  (V.-Vj) 
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The  name  Fourier  is  associated  with  the  first  term  on  the  right 


hand  side*  while  the  third  term  is  related  to  the  Dufour  effect.  The 

Q(lt^ 

depends  in  a  detailed  manner  on  the  frequency 
and  intensity  of  the  radiation  and  on  the  fluid  properties.  One 
procedure  that  may  be  used  at  present  is  to  evaluate  this  term  after 
having  solved  the  flow  field.  Obviously,  such  a  procedure  is  approximately 
valid  only  if  ^4  iQl  . 

The  work  done  by  the  pressure  tensor  is  given  by: 

'\[i  VV  Si  -  PV-V  +  $ 

-S.  I 


where  $  ,  the  dissipation  function  is  given  by: 


§  «  T;;  :  VV 


+ 


j  DXi  f 


\ 


y 


(31) 


(32) 


upon  having  introduced  equation  (24). 
For  example,  in  x,  y,  z  coordinates; 

”  ^ '  c)x  7>z  ' 


*  mt * [f,T *  (WT] 


(33) 
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Sattlag  X  ^  yU  and  raarraagiagi  the  diseipetioa  function  becomes: 


If  Stokes'  poetvlnte  holds,  $  %  O  •  That  is,  ^  is  essentially 
a  positive  quantity  diat  vanishes  only  ndien  the  fluid  motion  consists  of  a 
pure  idllatlon.  The  physical  meaning  Is  that  viscosity  is  a  method  of 
adding  Internal  energy  to  the  fluid  at  tiie  ejq  ense  of  other  forms  of  energy. 
For  a  noa«8tokeslan  fluid,  consisting  of  complex  molecules,  it  is 
possible  for  ^  <  O  . 

Introducing  equations  (30)  and  (31)  into  the  energy  equation,  the 
latter  becomes: 


pis  s  +  5  V*(K7T) 

dl'i 


4-  V 


.^21 

r\i  i  yi 


(35) 


The  energy  equation  may  also  be  written  in  terms  of  the  enthalpy. 

Taking 


equation  (33)  becomes: 
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or  expanded: 


=  v.(kvt^  -  V-Z  ^  ^  ^ 

dt  1  '  * 

+  v-EJ  22  iy,-V)  + 


If  we  consider  that  the  gaseous  system  is  composed  of  a  mixture  of  reacting 
perfect  gases«  then  the  caloric  equation  is  simply: 


■  (m 


c,  •  2  CiCvi 


Cl*  e,(T) 


e  •  ZciCi 


and  hence  upon  utilising  equation  (10),  the  energy  equation  may  be  written: 

/3Cy  dl  «  -7-Q  +  2/»iYi’^l 

d^  i 

-  2  eiWj  +  2  eiCv  piYi) 

which  may  be  expanded  into  the  form: 

f>Cv  ^  .  -  pv-v  +  §  +  V-{k^T)  -  V  'Zfiy^W^ 
dx  ^ 

*  '21^'  (V,-Vi)  +  SftYrfi 

^  *J  (4 

-  +  ?ei(v./5iVi) 

i  1  ' 
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We  can  also  obtain  a  similar  result  in  terms  of  the  enthalpy. 


We  take: 


Cfi  ■ 

* 

K(f) 

(42) 

Cp  • 

2  C-iCpi 

SCi^i 

and  there  follows  immediately: 

41  . 

dt 

v*(kvt)  - 

Yi-VT  +  $ 

+ 

V*  ?J 

njw  1 

V.-VA  ^  ^ 

dr 

nt  1 

1  -  J  -1/  (43) 

-ZhiWi  +  ZpiYi-fi 


It  should  be  noted  that  uhile  e  and  h  contain  the  chemical  er^r gy 
of  formatioat  tiiey  doa't  include  either  the  kinetic  energy  of  the  macroscopic 
medioa  of  Ihe  gas,  or  the  potential  energy  associated  with  the  external 
forces 

If  wa  udlh  to  include  the  kinetic  energy,  we  may  do  so  as  follows. 
Premultiply  the  momentum  equation  by  V  •  and  there  is  obtained: 


Oefins 

E  »  e  +• 

H  =  h-t-Av’-  »  e+P+.i^r’- 

/5  2 


(44) 


(45) 
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and  then  we  may  write  ^e  general  expreseione: 
# 

*  -v-Q  +  v-{Tr':jr)  -»• 


where 


v-dr  r)  =  ir 


Alto, 


*•  ^-(kvt)  -  v.ZpiYiK  +  ^ 

+  v  +  V.^ (V.  -va 

„  ^  n*  i),i 

+  f/’l'iiili 

We  have  obviously  not  exhausted  all  of  the  possibilities  for  representing 


the  energy  equation. 
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7.  THE  EQUATION  OF  STATE 


A  rtatoambla  approximation  to  the  thermal  behavior  of  the  gaseous 
mixturi’.  is  Sie  assumption  thig  the  eqnation  of  state  is  given  by: 


sAere 


p  *  ^  RT 

m 


(49) 


Thar*  alao  follows: 


R- 

A  • 


^  S  />iRi  - 

5T,  I 


2  CjRj 


2  Xi'Wi  * 


\ 


2  ^ 

i  ?ni 


(50) 


Soma  othar  uasful  form*  of  aquotioa  (49)  are: 

fi  • 

p  a 

Pt  •  NiAT 


Aflar  a  ■tafflels'At  aumbar  of  aolutlona  havs  boon  obtained  la  a 
glvaa  type  of  problorr.  utillaiag  par  fact  fa*  bshavlori  thsa  It  would  be 
wall  to  recoaaidar  the  effect  of  aon-ldeal  thermal  equation*. 
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8.  THE  INFLUSCID  EQUATIONS 

Up  to  thia  point,  all  of  the  eqiationt  have  been  rather  general. 

Let  ua  now  obtain  the  influacid  equations  by  neglecting  all  transport  phenomena. 
That  is,  diffusion,  viscosity  and  thermal  conductivity  are  eliminated 
from  further  consideration  in  the  inviscid  flow  field  and  there  follows 
accordingly; 

Conservation  of  species  i 


a 

<Li 

(52) 

Global  Continuity 

+  v  (pv’)  -  0 

H  ' 

(53) 

Momentum  Equation 

■  -vp  + 

(54) 

Energy  Equation 

pdK  -  dp 

pi^  »  ^  +  2  0  v  Fi 

at  st  1 

(55) 

which  shows  that  only  in  the  steady  state  is  the  stagnation  enthalpy 
.;on3tant  when  following  the  macroscopic  motioa  of  the  gas^  (neglecting 
external  forces). 
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In  tike  eolation  of  the  influecld  equntione,  farther  difficultiee  are 
introduced  -when  one  attempts  to  satisfy  the  boundary  conditions*  This 
situation  very  clearly  exists  in  a  hypersonic  flow  field  since  the 
governing  partial  differential  equations  are  of  mixed  type.  We  will  reserve 
our  discussion  until  later.  Let  us  now  proceed  to  the  boundary  layer 
equations. 
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9.  THE  BOUNDARY  LAYER  EQUATIONS 

The  first  thing  we  must  tske  cognisance  of  is  that  the  boundary 
layer  approximatioa  is  not  invariant  with  the  coordinate  system.  Stated 
another  way#  the  form  of  the  boundary  layer  equations  together  with  their 
boundary  conditions  actually  depend  on  the  choice  of  the  coordinate  system 
to  such  an  extent  that  different  results  may  actually  be  obtained  for  the 
same  physical  body  when  different  coordinate  systems  are  utilised^.  Thus» 
for  example«  if  the  boundary  layer  approximation  is  made  with  respect 
to  the  coordinate  system  sketched  in  figure  the  boundary  layer  equations 
for  a  reacting  multicomponent  mixture  can  be  formulated  as  follows « 

The  conservation  of  species  i  may  be  written: 

JP}  +  I  Pi(v =  Wi 
9t  -  -  - 

where  the  general  form  of  the  molar  diffusion  flux  is  given  by  eq.  (18) 

&1  °  ''lYl  ■  4  S  vJnT  (”) 

“  i'‘  iWi 

and  the  driving  force  ^  j  may  be  written: 

ij  «  vXj  +  (Xj  -  Cj)  vlnp  (58) 

If  •xtaroal  forces  are  neglected.  The  Dij  are  the  multicomponent 
diffusion  coefficients,  of  which  there  are  where  M  is  the 

number  of  species,  and  the  X 


are  the  mole  fractions. 
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Defining  a  referred  mole  fraction  Bi « 
diifttsion  equation  may  be  written: 

^  +  V-  j),  -  M. 

dt  4-* 

On  introd  Xing  eqe«  (57)  and  (58),  eq.  (59)  becomes: 


the 


(59) 


When  the  boundary  layer  approximation  it  made  with  retpect  to  a 
body-oriented  axially-tymmetric  coordinate  system,  the  diffusion  equation 


becomes: 

O 

e>X 

+  u'5li  2!i  Lii5.(£BiV!dir)1.M  (6i) 


pirn" 


The  energy  equation  may  be  written: 


a  ^  +  v{<VT)  -  2  Cp  /OVi-VT-  Z  loiVV^  (62) 

dt  i  »  ^ 

where  the  Dufour  effect  has  been  neglected  and  the  following  notation  has 
been  introduced. 


Vii*  aK/,  +  fJcp.JT;  c^-ZciCfi  («) 

On  making  the  boundary  layer  approximation  there  is  obtained: 

-  Mj  I-,  (*»j)  ■  i;a  i  a  -  2 


I 
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It  is  again  noted  that  equation  (62)  merely  represents  one  of  several 
convenient  forms. since  the  energy  equation  can  be  written  in  many 
equivalent  forms*  Global  continuity  becomes: 


Finally  upon  neglecting  external  forces  the  momentum  equation  yields  the 
following  X  and  y  scalar  equations: 


The  behavior  of  the  above  equations^  under  the  influerxe  of  suitable  I 

constraints,  forms  the  basis  for  studies  of  the  interaction  between  a  j 

surface  and  a  reacting  stream* 

It  is  noted  that  under  certain  circumstances,  it  may  be  permissible 
to  lump  the  gaseous  components  whose  thermal  and  chemical  behavior  is 
similar  so  that  although  the  gas  is  actually  a  multicomponent  mixture,  its 
behavior  may  be  approximated  by  a  binary  mixture.  Further,  in  formulating 
the  governing  equations  for  a  two -component  chemically  reacting  boundary 
layer,  it  may  be  convenient  to  utilize  mass  fractions  rather  than  mole 
fractions  in  representing  the  composition  of  the  gas. 

I 

1 

For  a  binary  mixture  wc  can  introduce  .q,  (Zi)  so  that  the  diffusion  i 

I 

i 

i 
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and  energy  equations  become  respectively: 


al  t  u-sci  +v3ti\ 
'  ax  au  ' 


(68) 


while  the  global  continuity  and  momentum  equations  remain  unchanged.  Here 


agaiOg  all  possible  forms  of  the  energy  equation  have  not  been  exhausted. 


For  esuunpleg  another  useful  form  is  obtained  from  eq.  (38) 


Note  that  the  Dufour  effect  and  all  external  forces  have  been  takers 
negligible  above. 

Ws  have  also  assumed  that  we  have  a  sufficient  number  of  auxiliary 
relations  for  the  determination  of  the  transport  properties  of  the  gaseous 
mixture.  For  example,  to  a  good  approximation  the  viscosity  of  a  mixture 
in  terms  of  the  viscosities  of  the  pure  components,  the  binary  diffuiion 
coefficients  and  the  mole  fractions  is  given  by*: 


1.385 


r- 

P  .OSk 


(«) 
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10.  BOUNDARY  CONDITIONS 

The  unknown  dependent  variables  in  a  reacting  flow  field  are  the 
three  velocity  components  t  (this  reduces  to  two  for  axially  symmetric 
flow)t  the  pressure^  the  temperature  and  the  composition.  Once  these  have 
been  determined^  the  problem  may  be  considered  solved.  ThuSt  in 
generalt  the  mathematical  problem  consists  of  solving  for  N  5  unknowns^ 
where  N  is  the  total  number  of  chemical  species.  To  accomplish  this 
endg  we  have  at  our  disposal  the  following  system  of  equations  ?  N  4  5 
in  total! 

N  -  1  conservation  of  species 
1  global  continuity 

3  momentum  equations 

1  energy 

1  state 

In  additiong  we  must  also  provide  a  set  of  boundary  conditions 
which  is  in  agreement  with  the  overall  order  of  the  mathetnatical  system. 
Examination  shows  that  the  order  of  the  system  ani  hence  its  mathematical 
character  depends  on  whether  the  flow  field  is  fluscid  or  influscid.  In 
generalg  when  the  flow  field  is  influscid,  the  system  is  of  the  N  +  4^  order, 
however,  when  the  flow  is  fluscid,  the  overall  order  of  the  system  is  2N  +  7, 
and  hence  the  fluscid  system  requires  N  +  3  more  boundary  conditions 
than  the  influscid  flow  system,  in  order  to  properly  define  the 
mathematical  problem. 
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The  relevant  boundary  conditions  or  compatibility  relations  for 
a  reacting  boundary  layer  must  be  developed  individually  for  each  problem* 

It  is  remarked  that  the  difficulty  at  the  surface  usually  lies  in  the  coupled 
mass«energy  transport*  Typical  treatments  of  the  surface  boundary  conditions 
for  a  reacting  boundary  layer  appear  in  references  4  and  5*  The  boundary 
conditions  at  the  outer  edge  of  the  boundary  layer  must  be  obtained  from 
the  influscid  solution*  It  is  therefore  of  the  greatest  importance  that 
solutions  be  obtained  for  the  non«equilibrium  influscid  flow  fields  since 
unless  one  knows  the  conditions  at  the  outer  edge  of  the  boundary  layer » 
oat  eaanet  prsdlet  the  state  of  affairs  inside  the  fluscid  layer. 
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11.  CONCLUSIONS 

Our  major  interest  here  has  been  in  formulating  the  boundary 
layer  equations  and  not  in  determining  heat  and  mass  transfer  which  will 
be  considered  elsewhere* 

From^^^  brief  discussion,  it  is  clear  that  the  analytical  treatment 
of  flow^e^s  in  which  chemical  reactions  are  present  requires  a  major 
e^ljli^with  regard  to  the  solution  of  the  governing  equations,  subject  to  the 
/appropriate  boundary  conditions. 


y 


While  general  techniques  for  obtaining  closed  form  solutions  of 
the  partial  differential  equations  are  desirable,  the  advent  of  high  speed 
digital  and  analog  computers  has  'made  it  more  likely  that  particular 
numerical  solutions  will  be  obtained,  and  hence  it  is  anticipated  that 
experiments  will  have  to  be  carefully  designed  to  detect  the  effects  which 
can  be  treated  theoretically. 
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COORDINATE  SYSTEM  FOR  AN  AXIALLY  SYMMETRIC  BODY 


Figure  I 


SCHEMATIC  OF  STAGNATION  POINT  BOUNDARY 

LAYER 


Figure  3 


